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Abstract 

Under the assumption that the dijet excess seen by the CDF Collaboration near 
150 GeV in Wjj production is due to the lightest technipion of the low-scale technicolor 
process pt — ^ Witt, we study its observability in LHC detectors for y/s = 8TeV 
and J Cdt — 20 fb -1 . We describe interesting new kinematic tests that can provide 
independent confirmation of this LSTC hypothesis. We show that cuts similar to those 
employed by CDF, and recently by ATLAS, cannot confirm the dijet signal. We propose 
cuts tailored to the LSTC hypothesis and its backgrounds at the LHC that may reveal 
Pt — > £vjj- Observation of the isospin-related channel p^ —> Ztt^ — > and of 

Pt — > WZ in the l+^l^vn and modes will be important confirmations of the 

LSTC interpretation of the CDF signal. The Ztvt channel is experimentally cleaner 
than Wtvt and its rate is known from Wtvt by phase space. It can be discovered or 
excluded with the collider data expected by the end of 2012. The WZ — > 3£v channel 
is cleanest of all and its rate is determined from Wttt and the LSTC parameter sinx- 
This channel and WZ —> £ + £~jj are discussed as a function of sinx- 
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1. Introduction 



The CDF Collaboration has reported evidence for a resonance near 150 GeV in the dijet- 
mass spectrum, Mjj, of Wjj production. This was based on an integrated luminosity of 
4.3 fb -1 pQ and updated with a total data sample of 7.3 fb -1 [2]. In Ref. [2], the resonant 
dijet excess has a significance of 4.1 a. The D0 Collaboration, on the other hand, published 
a search for this resonance based on 4.3 fb -1 that found no significant excess. Based on a 
ty+Higgs boson production model, D0 determined a cross section for a potential signal of 
O.82I0I2 Pb an d a 95% confidence level upper limit of 1.9 pb [3]. Analyzing its data with the 
same production model, CDF reported a signal rate of 3.0±0.7pb and a discrepancy between 
the two experiments of 2.5 a [4J. This discrepancy remains. The purpose of this paper is 
to help guide the LHC experiments in searches to test for the CDF dijet excess in the Wjj 
and two closely related channels. We do this in the context low-scale technicolor (LSTC), 
interpreting CDF's dijet excess as the lightest technipion np of this scenario, produced in 
association with W ± in the decay ap° —> Wttt and decaying to a pair of quark jets j5]. 
The related channels supporting this interpretation are p^,a^ —> Zn^ and W^Zj^] They 
require no additional LSTC model assumptions beyond those made in Ref. [5j to determine 
LHC production rates. We assume y/s = 8TeV and consider J Cdt = 20 fb -1 , the amount 
of data expected to be in hand by the end of 20120 

Low-scale technicolor (LSTC) is a phenomenology based on walking technicolor 0, El 
[TUl [TT] . The gauge coupling qitc must run very slowly for 100s of TeV above the TC scale, 
A T c ~ several 100 GeV, so that extended technicolor (ETC) can generate sizable quark and 
lepton masses while suppressing flavor-changing neutral current interactions [12] . This may 
be achieved, e.g., with technifermions belonging to higher-dimensional representations of the 
TC gauge group. Then, the constraints of Ref. [12] on the number of ETC-fermion repre- 
sentations imply that there will be technifermions in the fundamental TC representation as 
well. They are expected to condense at an appreciably lower energy scale than those belong- 
ing to the higher-dimensional representations and, thus, their technipions' decay constant 
F\ « = (246 GeV) 2 [T3] . Spin-one bound states of these technifermions will have an 
orthoquarkonium-like spectrum with masses well below a TeV — greater than the previous 
Tevatron limit M PT ^ 250 GeV [THUS] and probably less than 600-700 GeV, a scale at which 
we believe the notion of "low-scale" TC ceases to make sense. The most accessible states are 
the lightest techni vectors, Vt = Pt{I G J PC = 1 + 1 ), Ut(0~1 ) and aT(l~l ++ ). Through 
their mixing with the electroweak bosons, they are readily produced as s-channel resonances 
via the Drell-Yan process in colliders. Spin-zero technipions ttt(1~0~ + ) are accessed in Vt 

1 LHC studies of the Wjj and WZ channels carried out so far are discussed in Sees. 3 and 5, respectively. 

2 Preliminary versions of this paper were circulated in Ref. [6 and Ref. [7 assuming y/s = 7 TeV and 
J Cdt = 1-20 fb _1 . The simulations in the current paper may be applied to different luminosities by scaling 
the event rates. We have not included the nontrivial effects of pileup at the higher luminosities of 8-TeV 
running. They also make difficult a detailed comparison of our results with the earlier 7- TeV ones. Our 
signal cross sections are uniformly 20% greater at 8 TeV than at 7 TeV, but the increases in various physics 
backgrounds are not so simply summarized. 
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decays. A central assumption of LSTC is that these lightest technihadrons may be treated 
in isolation, without significant mixing or other interference from higher-mass states. Also, 
we expect that (1) the lightest technifermions are 5C/(3)-color singlets, (2) isospin violation 
is small for V T and 7Tt, (3) M UT = M PT) and (4) M aT is not far above M PT . This last 
assumption is made to keep the low-scale TC contribution to the S^-parameter small. An 
extensive discussion of LSTC, including these points and precision electroweak constraints, 
is given in Ref. [T6] . 

Walking technicolor has another important consequence: it enhances M nT relative to 
M PT so that the all-7Tr decay channels of the Vt are likely to be closed [13J. Principal Vp- 
decay modes are Witt, Zttt, ^ttt, a pair of EW bosons (which can include one photon), 
and fermion-antifermion pairs [T71 HHJ [[6]. If allowed by isospin, parity and angular mo- 
mentum, Vt decays to one or more weak bosons involve longitudinally-polarized Wl/Zl, 
the technipions absorbed via the Higgs mechanism. The rates for these nominally strong 
decays are suppressed by powers of sin 2 \ = (Fi/F^) 2 <C 1. This important LSTC parameter 
is a mixing factor that measures the amount that the lowest-scale technipion is the mass 
eigenstate ttt (cosx) and the amount that it is Wl/Z l (sin%). Thus, each replacement of 
a mass-eigenstate ttt by Wl/Z l in a 1^ decay amplitude costs a factor of tan^- Decays 
to transversely-polarized 7, W±,Zj_ are suppressed by g,g'. Thus, the Vp are very narrow, 
T(pt) ^ 1 GeV and T^t, clt) ^0.1 GeV for the masses considered here. These decays have 
striking signatures, visible above backgrounds within a limited mass range at the Tevatron 
and probably up to 600-700 GeV at the LHC [111 [20]. 

In Ref. [5j we proposed that CDF's dijet excess is due to resonant production of Wttt 
with M nT = 160 GeV. We took M PT = 290 GeV and M aT = 1AM PT = 320GeVQ Then, 
about 75% of the Wttt rate at the Tevatron is due to Pt —> Wttt and, of this, most of the 
Ws are longitudinally polarized^] The remainder is dominated by production. Its decay, 
and a small fraction of the Pt's, involve W± production, which is generated by dimension-five 
operators [16] . These operators are suppressed by mass parameters M v ,a that we take equal 
to M PT . The other LSTC parameters relevant to Wttt production are g PT7TT7TT and sin;^- 
The pt — >• tttttt coupling g PT7TT7TT is the same for all Pt decays considered here and it is 
naively scaled from QCD; its Pythia default value is gl p rp — g * n /4tt = 2.16(3/A^ TC ) with 
N T c = 4. We use sinx = 1/3. Using the LSTC model implemented in Pythia fT71 [T8l I2T] . 
we found a(pp -> p T -> Wttt Wjj) = 2.2 pb (480 fb for W ev, Adopting CDF's 

cuts, we closely matched its Mjj distribution for signal and background. Motivated by the 
peculiar kinematics of Pt production at the Tevatron and Pt — >• Wttt decay, we also suggested 



3 The Pythia default decays for technipions are based on the assumption that they are Higgs-like, i.e., 
involve couplings proportional to fermion mass. They are thus dominated by 7rJ — » c6, ub and tt§, — >> bb. 
These modes involve energy loss to neutrinos that we have not included in reconstructing dijet masses. 
Therefore, the choice M nT = 160 GeV reconstructs close to 150 GeV. If technipions decay mainly to light 
quarks and leptons, a plausible possibility for the lightest 7Tt, then we would expect all our input technihadron 
masses to decrease by 10-15 GeV. 

4 About 70% of the Wttt rate at the LHC is due to the pr- 

5 This includes B(tt t -> qq) ^ 90% in the default Pythia 7r T -decay table. 
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cuts intended to enhance the ttt signal's significance and to make Pt —> Wjj visible. Several 
distributions of data in the excess region 115 GeV < Mjj < 175 GeV published by CDF [2] 
— notably Mwjj, Pr(jj), A0 an d Ai? — \/ (Arj) 2 + (A0) 2 — fit the expectations of the 
LSTC model very well. The background-subtracted AR distribution, in particular, has a 
behavior which, we believe, furnishes strong support for our dijet production mechanism. 

The purpose of this paper is to propose and study ways to test for the CDF signal at 
the LHC. In Sec. 2 we review the kinematics of Pt^t — >• Wttt and Zn T in LSTC. We also 
present an interesting new result: the nonanalytic behavior of da/d(AR) and da/d(Ax) 
at their thresholds, (Ai?) min and (Ax)min- Here Ax is the opening angle between the ttt 
decay jets in the Pt rest frame. For massless jets, a good approximation, we find that 
(Ai?) min = (Ax)min = 2cos _1 (^), where v = r p^ T jE^ T is the ttt velocity in the Pt rest 
frame. This result, peculiar to production models such as LSTC in which a narrow resonance 
decays to another narrow resonance plus a W or Z, provides measures of v independent of 
p/E and, hence, valuable corroboration of this type of production. In Sec. 3 we consider 
the Pt,o>t Wttt process. Its LHC cross section at 8TeV is 9.5 pb but, for CDF cuts, 
its backgrounds have increased by about a factor of ten over those at the Tevatron. This 
makes testing for the dijet excess in this channel very challenging. We suggest cuts which 
enhance signal-to-background (S/B) but which will still require a very good understanding 
of the backgrounds in Wjj production. Recent studies of Wjj production by ATLAS and 
CMS are discussed there. In Sec. 4 we study p^,a^ — >> Ztt^ whose cross section is 2.8 pb 
at 8TeV (190 fb after Z —> e + e~, p + p~). This is the isospin partner of p^,a^ Wir^, so 
its cross section is rather confidently known. The £ + £~jj channel is free of QCD multijet 
and it backgrounds and missing energy uncertainty. Reconstructing the Zjj invariant mass 
and other signal distributions, particularly in AR and Ax, will benefit from this. Because 
of these features, we believe that the Zn T — >> Zjj mode will be the surest test of CDF's dijet 
signal at the LHC. In Sec. 5, we study p^,a^ — >> WZ. The cross section for this mode is 
proportional to tan 2 x times the p^ , a? ~~ ^ W ±/ k { ^ and Ztt^ rates, but enhanced by its greater 
phase space. We predict a(p^,a^ — >> WZ) = 1.8 (1.1) pb for sinx = 1/3(1/4). In the all- 
leptons 3£v mode with e's and /x's, the rate is only 26(15)fb, but jet-related uncertainties 
are absent except insofar as they effect $t resolution. A new study by CMS of this channel 
is discussed there. The WZ — >> £^£~jj mode is also an interesting target of opportunity so 
long as sinx ^ 1/4. The AR and Ax distributions for Z —> jj again provide support for 
our narrow LSTC-resonance production model. In short, one or both of the Zn T and WZ 
modes should be dispositive of the LSTC interpretation of the CDF dijet excess with the 
~ 20 fb -1 expected by the end of 2012. We present in an appendix the details of calculations 
in Sec. 2 regarding the nonanalytic threshold behavior of the Ax and AR distributions. 

While the simulations of the CDF signal in this paper are made in the context of low- 
scale technicolor, their qualitative features apply to any model in which that signal is due 
to qq production of a narrow resonance decaying to a W plus another narrow resonance. 
Several papers have appeared proposing such an s-channel mechanism [22l [23j EH [25l [26j [27] . 
With similar resonance masses to our LSTC proposal, these models will have kinematic 
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distributions like those we describe in Sec. 2. However, not all these models will have the 
Zjj and WZ signals of LSTC. There are also a large number of papers proposing that 
the CDF signal is due to production of a new particle (e.g., a leptophobic Z') that is not 
resonantly produced [2HI [29l [30l EH [32l [33l E3] . These "t-channel" models will not pass our 
kinematic tests. 

2. LSTC Kinematics and Threshold Nonanalyticity 

The kinematics of Pt — >• Wttt at the Tevatron and LHC are a consequence of the basic LSTC 
feature that walking TC enhancements of M^ T strongly suggest M PT < 2M^ T and, indeed, 
that the phase space for p T Wttt is quite limited [13l [35] . At the Tevatron, a 290 GeV pr 
is produced almost at rest, with almost no Pt and very little boost along the beam direction. 
At the LHC, Pt(pt) ^ 25 GeV and 77 (py) ^ 2.0. Furthermore, the ttt is emitted very slowly 
in the Pt rest frame — v ~ 0.4 for our assumed masses — so that its decay jets are roughly 
back-to-back in the lab frame. Thus, Pt(ttt) & 80 GeV and the z-boost invariant quantities 
Acj) and AR = y/ (Arj) 2 + (Acj)) 2 are peaked at large values less than tt. 

These features of LSTC are supported by CDF's 7.3 fb" 1 data [2j. Fi gures 1H4 show 
distributions before and after background subtraction taken from the 115 < Mjj < 175 GeV 
region containing the dijet excess. The subtracted-data M W jj signal has a narrow resonant 
shape quite near 290 GeV. Unfortunately, the background peaks not far below that mass so 
that one may be concerned that the subtracted data's peak is due to underestimating the 
background. Also, as we expect, the subtracted pr(jj) data falls off sharply above 75 GeV 
and the subtracted Acj) data is strongly peaked at large values. Again, one may worry that 
these are artifacts of the peak of the Mwjj background and the position of the Mjj excess. 

The background-subtracted AR distribution, however, is very interesting. It is practically 
zero for AR < 2.25, then rises sharply to a broad maximum before falling to zero again at 
AR ~ 3.5. This behavior, and a somewhat similar one we predict for Ax are the main subject 
of this section. We will show that the threshold form of the AR and Ax distributions provide 
direct measures of the velocity of the dijet system in the subprocess center-of-mass frame that 
are independent of measuring p/E and, thus, are independent checks on the two-resonance 
topology of the dijet 's production mechanising One might think that the corresponding 
ARa and A\u distributions from Z —> £^£~ would be similarly valuable. Unfortunately, 
because the dileptons come from real Z's and our cuts make the background Z's like the 
signal ones, ARu and A\u are indistinguishable from their backgrounds. 

For our analysis, we assume the jets from tit decay are massless. We have examined the 
effect of including jet masses and found them to be unimportant. We will remark briefly on 
this at the end of this section. We first consider the dominant Pt contribution to W/Zttt 

6 Note that AR and A% are largely unaffected by lost neutrinos if semileptonic 6-decays are an important 
component of ttt decays. Also, A% is defined in the pr rest frame, while AR is defined in the lab frame. If 
one wishes to remove the effect of Pt(Pt) on Ai?, it should be defined in the pT frame. 



5 



CDF Run II Preliminary L int = 7.30 fb 1 



CDF Run II Preliminary L = 7.30 fb' 1 

J int 



0> 



LU 



600- 



400- 



200- 



0, 




200 400 600 

KS = 73.5 %, % /ndf= 1 44.5/ 46 M QeV/c 2 

Ivj] 



03 

Q 



100 




100 200 300 400 500 600 

mtot 



Figure 1: CDF Mwjj distributions for J Cdt = 7.3 fb -1 from the dijet signal region 115 < 
Mjj < 175 GeV J2]. Left: Expected backgrounds and data; right: background subtracted 
data. 
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Figure 2: CDF pr(jj) distributions for J Cdt = 7.3 fb -1 from the dijet signal region 115 < 
Mjj < 175 GeV [2]. Left: Expected backgrounds and data; right: background subtracted 
data. 
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Figure 3: CDF Acj) distributions for J Cdt = 7.3 fb -1 from the dijet signal region 115 < 
Mjj < 175 GeV [2]. Left: Expected backgrounds and data; right: background subtracted 
data. 
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Figure 4: CDF AR distributions for J Cdt = 7.3 fb -1 from the dijet signal region 115 < 
Mjj < 175 GeV [2]. Left: Expected backgrounds and data; right: background subtracted 
data. 
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production, commenting on the contribution also at the end. 

Define the angles 0, #* and 0* as follows: Choose the z-axis as the direction of the event's 
boost; this is usually the direction of the incoming quark in the subprocess cm. frame. 
In the Pt rest frame, 9 is the polar angle of the ttt velocity the angle it makes with 
the z-axis. Define the xz-plane as the one containing the unit vectors z and v ) so that 
v = x sin 9 + z cos 9 ) and y — zxx. Define a starred coordinate system in the ttt rest frame 
by making a rotation by angle 9 about the ?/-axis of the Pt frame. This rotation takes z 
into z* = v and x into x* = xcos9 — zsin#. In this frame, let p[ be the unit vector in the 
direction one of the jets (partons). The angle between v and p* is the azimuthal angle 
of p* = -p* is 0*: 

cos # = cos9*=pl'V, tan (/>* = pl y * / p[ x * . (1) 

Note that, since ttt —t qq is isotropic in its rest frame, da(qq — >> Pt —> W j j) / d(cos 9*) = a/2, 
where a is the total subprocess cross section. 

It is easier to consider the da/d(Ax) distribution first. For massless jets, 

1 - cos(Ax) = 1 2(1 2 ~^ • (2) 
1 — v z cos 2 9* 

The minimum value of Ax occurs when #* = n/2 (i.e., v _L p*), and so 

7T > A X > ( Ax) min = 2 cos" 1 (v) . (3) 
From Eq. ^ , it is easy to see that 

da (1 — v 2 ) a 

(4) 



d(A X ) 4^sin 2 (Ax/2) Vcos 2 ((Ax) min /2) - cos 2 ((Ax)/2) ' 

The Ax distribution has an inverse-square-root singularity at Ax = (Ax) m in = 2cos _1 (^) = 
2.23 for our input masses, and falls sharply above there. This is illustrated in Fig. [5] where 
we plot this distribution for the primary partons and for the reconstructed jets. The low-side 
tail for the jets is an artifact of their reconstruction. 

To understand this singularity better, it follows from Eq. ^ that Ax may be expanded 
about cos 0* = as 

Ax= (Ax) min + ^cos 2 r +-.. , (5) 

where a is a positive ^-dependent coefficient. Then, near cos#* = 0, i.e., the Ax threshold, 

da a d(cos9*) 1 , N 

~ <x , = • (6) 



rf(Ax) 2 rf(Ax) ^Ax - (Ax) 

It is the simple one-variable Taylor expansion of Ax in Eq. ^ that has caused this singu- 
larity. 

The discussion of da/d{AR) for the LSTC signal shares some features with da/d{A\) 
though it it is qualitatively different. The AR distribution also vanishes below a threshold 
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Ax in p T rest frame, area normalized AR- in p T rest frame, area normalized 




Figure 5: The area-normalized Ax and AR distributions for the primary parton/jet in 
Pt,Q>t - ^ Wttt production followed by ttt —> qq decay, constructed as described in the text. 
Red: pure distribution of primary parton before any radiation; blue: the distribution for the 
jets reconstructed as described in Sec. 3. 



(Ai2) m i n , which is equal to (Ax)min = 2cos _1 ('u). This remarkable feature, derived in the 
appendix, can be understood simply as a consequence of the fact that the minimum of AR 
occurs when both jet rapidities vanish. In that case, AR = A(j) = A%. 

At threshold, however, the AR distribution is oc ^ AR — (Ax)min, not the inverse square 
root. As illustrated in Fig. [5| it rises sharply from threshold into a broad feature before 
decreasing. The measure of the ttt velocity v is given by the onset of the rise, not its 
peak. This is the behavior seen in the CDF data in Fig. [4j where the rise starts very near 
2cos _1 ('u) = 2.23 for our input masses. Both the Ax and AR distributions measure the 
ttt velocity v and, therefore, provide confirmations of the Pt — >• Wttt hypothesis which are 
independent of the background under the M W jj resonant peak and of uncertainty in the $t 
resolution as well. 

The reason for this qualitative difference between the two distributions is that da/ Ax 
involves a one-dimensional trade of cos#* for Ax, whereas AR is parametrized in terms 
of the three angles 6,6*, 0* in an intricate way, with all three being integrated over to 
account for the constraint defining AR. In contrast to what happens in the Ax case, the 
Jacobian singularity at the threshold is "antidifferentiated" twice, hence its comparatively 
lower strength. Using a Fadeev-Popov-like trick, the AR distribution can be written 

rf(cos 6) rf(cos 0*)d(cos 0*) - da - r 5 (AR - /(cos 6, cos6*, cos (/>*)) . (7) 
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The function /(cos 0, cos 0*, cos 0*) is shown in the appendix to have its absolute minimum 
at cos# = cos#* = cos0* = 0, for which its value is equal to (Ax)min- Near its minimum it 
is locally parabolic and its Taylor expansion is 

/(cos 0, cos 0*, cos 0*) = (Ax)min + \ (bo cos 2 9 + b e * cos 2 0* + b r cos 2 (jf) H (8) 



The positive ^-dependent coefficients bo, bo* and b^ are also given in the appendix, Eq. (29). 
For AR close to (Ax) nun, this expansion can be used to approximate Eq. ([7|. In a similar way 
as for the Ax distribution, integrating first over cos 9* generates the appearance of a Jacobian 
inverse square root singularity oc [2(AR — (Ax)min) ~~ (bo cos 2 9 + b^ cos 2 0*)] -1 / 2 . The two 
remaining integrations over cos 9 and cos 0* were trivial in the Ax case as the integrand did 
not depend on them, but this is not so for AR which involves a double integration over a 
restricted angular phase space defined by 

< b cos 2 9 + V cos 2 0* < 2 (AR - (Ax) min ) • (9) 

Performing the integral in Eq. ^ near (Ai?) min = (Ax) m in yields a result oc AR — (Ax)min- 
We have examined the effect of finite jet masses (as opposed to jet reconstruction and 
energy resolution) on the threshold values of the AR and Ax distributions and the extraction 
of the ttt velocity v from them. Our jets (which include 6-jets in the Pythia default 7Pr-decay 
table) have masses ^ lOGeV. Assuming, for simplicity, equal jet masses and denoting by 
u = yjl — 4Mj 2 t /M^ T the jet velocity in the ttt rest frame, the corrected (Ax)min(^) is 

(AxWu) = cos" 1 v2 ~ U l i { \~ V 2 - cos" 1 ^ 2 - 1) - v(l - v 2 yi 2 (l - u 2 ) . (10) 
v z + ir (1 — v z ) 

This is less than the massless (Ax) m in by half a percent for Mj et = lOGeV. 

Finally, as noted, the accounts for about 25-30% of Wttt production. This decay 
gives a ttt velocity of 0.54 in the rest frame and (Ax) m in = 2.00. The effect is clearly 
visible in the Ax and AR distributions for the primary parton in Fig. [5j but is washed out 
by the low-end tails for the reconstructed jets. We believe that the low and high-end tails 
are due to the two ttt jets fragmenting to three jets and the two leading jets being closer 
or farther apart than the original pair. It turns out that our Q- value cut for Zttt in Sec. 4 
eliminates the contribution to the signal. 



3. The p T , a T — ^ Wtt t mode at the LHC 

As a reminder, we assumed M PT = 290 GeV, M aT = 1.1M PT = 320 GeV, M nT = 160 GeV 
and sinx = 1/3 to describe the CDF dijet excess. The Tevatron cross section is 2.2 pb. At 
the 8 TeV LHC, these parameters give g{Wtit) = 9.5 pb (2.0 pb for W —> ev, \iv). These 
cross sections are 20% higher than at 7 TeV, but this does not translate into a 20% increase 
in S/B. About 70% of the LHC rate is due to the pr] the pr and ar interference is very 
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Figure 6: Left: The ATLAS Mjj distribution for exactly two jets in Wjj production at 
y/s = 7TeV and J Cdt = 1.02 fb -1 ; from Ref. [36] Right: Simulation of the Mjj distribution 
in Wjj production with ATLAS cuts (except that Pt{^) > 30GeV) for 1.0 fb -1 . The open 
red histogram is the ttt —> jj signal times 10. 



small. For such close masses, it is impossible to resolve the two resonances in the Mwjj 
spectrum. 

Last summer, the ATLAS Collaboration published dijet spectra for 1.02 fb -1 of Wjj data 
with exactly two jets and with two or more jets passing selection criteria [36]. The ATLAS 
cuts, taken as close to CDF's as practical, were: one isolated electron with E T > 25GeV or 
muon with p T > 20GeV and rapidity \r] £ \ < 2.5; $ T > 25GeV and M T {W) > 40GeV; two 
(or more) jets with p T > 30GeV and \rjj\ < 2.8; and Pr(jj) > 40GeV and Ar] < 2.5 for 
the two leading jets. The Mjj distribution for the two-jet data is shown in Fig. [HJ There is 
no evidence of CDF's dijet excess near 150 GeV nor even of the standard model WW/WZ 
signal near 80 GeV. This is what we anticipated in Ref. [6] because of the great increase in 
Wjj backgrounds at the LHC relative to the Tevatron. On the other hand, it is noteworthy 
and encouraging for future prospects that the ATLAS background simulation appears to fit 
the data well. 

In Fig. [6] we also show our simulation of the LSTC Mjj signal and backgrounds at 
the LHC for = 7TeV and J Cdt = 1.0 fb -1 . ATLAS's cuts were used except that we 
required pr(^) > 30GeVj^] This tighter cut and our inability to include the data-driven 

backgrounds were generated at matrix-element level using ALPGENv213 [37] . then passed to 
Pythiav6.4 for showering and hadronization. We use CTEQ6L1 parton distribution functions and a factor- 
ization/renormalization scale of \i = 2Mw throughout. For the dominant W+jets background we generate 
W + 2j (exclusive) plus W + Sj (inclusive) samples, matched using the MLM procedure [38 (parton level 
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QCD background account for our lower event rate compared to ATLAS. Despite this, the 
agreement between the two is quite good. In particular, our simulation shows that the 
CDF/ ATLAS cuts can neither reveal nor exclude the LSTC interpretation of the CDF signal 
at the LHC for any reasonable luminosity^] 

Recently, the CMS Collaboration studied the dijet-mass spectrum in W(—t tv) plus 
jets production with 4.7 fb -1 at 7TeV [43]. CMS used the following cuts which were partly 
adopted from Ref. [7J: p T (e, fj) > 25,30GeV and rapidity |??(e,/x)| < 2.5,2.1, AR(£J) > 0.3; 
#r(e,/x) > 35,25GeV, A</>($ T J) > 0.4; M T (W) > 50GeV and p T (W) > 60GeV; exactly 
two or three jets with p T i > 40GeV, p T 2,3 > 30GeV, \rjj\ < 2.4; and Pr{jj) > 45GeV, 
Ar)(jj) < 1.2. CMS used MadGraph to generate W + jets and a data-driven method 
to determine the Mjj shape and background: A superposition of a set of templates was 
constructed in which the MadGraph factorization and renormalization scales were varied 
up and down by a factor of two from their default values, and this was fit to the dijet 
spectrum outside the signal region, taken to be 123 to 186 GeV. The Wjj background in 
the signal region was then determined from this fit. The CMS dijet spectra before and after 
background subtraction are shown in Fig. [7| Note that the vertical scale is a Events/GeV." 
No significant enhancement near 150 GeV was observed. (What CMS meant by a a CDF-like 
signal" is not specified in Ref. [43].) Using a WH production model, CMS reported a 95% 
upper limit on the production cross section times B(W — >• tv) of 1.3 pb. 

We studied the LSTC Wjj signal at yfs = 7TeV in Ref. [7J, before the CMS paper's 
release. Our prediction for the cross section was aB = 1.7 pb, 30% higher than CMS's limit. 
In order to achieve a better outcome than ATLAS 's 2011 study, we examined a variety of 
cuts motivated by Pt —> Wttt kinematics. Cuts quite similar to those we proposed for the 
Tevatron in Ref. [5] typically caused the background to peak very near the dijet resonance. 
To get the signal off the peak (and more like the original CDF Mjj excess [I]), we used 
the following: lepton p Tl > 30 GeV and \rj £ \ < 2.5, $ T > 25 GeV, M T (W) > 40 GeV 
and p T (W) > 60 GeV; exactly two jets with p T1 > 40 GeV, p T2 > 30 GeV, \rjj\ < 2.8; 
Pr(jj) > 45 GeV, Arj(jj) < 1.2; and Q = M Wjj - M i5 - M w < 100 GeV. The resulting M j5 
distribution is also displayed in Fig. [7j Counting events in the range 120 < Mjj < 170 GeV 
gives S/y/B = 6.5 for this luminosity, but only S/B = 0.050. The AR and Ax signals are 
also small and not useful. Because of the small S/B, and in view of the difficulty CMS had 
fitting the dijet spectrum in the diboson and CDF-signal region, we believe that a better 



cuts are imposed to ensure that W + 0, 1 jet events cannot contribute). After matching, the overall normal- 
ization is scaled to the NLO W + jj value, calculated with MCFMv6 [39]. After passing through Pythia, 
final state particles are combined into (77, 0) cells of size 0.1 x 0.1, and the energy in each cell smeared with 
AE/E = 1.0/y/E/GeV. The energy of each cell is rescaled to make it massless. Isolated photons and 
leptons (e, /i) are removed, and all remaining cells with energy greater than 1 GeV are clustered into jets 
using Fast Jet (anti-kT algorithm, R = 0.4) [40]. Estimates of the background including higher order effects 
have been shown to be completely consistent with our LO+PS treatment [4TJ [42] . Finally, the quadratic 
ambiguity in the W reconstruction was resolved by choosing the solution with the smaller p z (v). 

8 Models of the CDF signal that are gg-initiated or involve large coupling to heavier quarks, e.g., Refs. [27j 
132] . are likely excluded by the ATLAS data. 
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CMS preliminary, L = 4.7 fb -1 , y/s = 7TeV 
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Figure 7: The CMS Mjj distributions for 4.7 fb -1 of W \lv, ev plus two or three jets 
data at yfs = 7TeV before (top left) and after (top right) the background subtraction 
summarized in the text; from Ref. [43J. On the bottom is our Mjj distribution for the 
Ptj a r —> Wttt —> £vdj signal and backgrounds at the LHC for 5fb _1 . Augmented ATLAS- 
like cuts as described in the text were used. The open red histograms are the ttt and Pt 
signals times 10. 



understanding of the backgrounds is required to observe or exclude the LSTC signal in this 
channel. 

Our simulations of the Mjj and M W jj distributions in Wjj production at y/s = 8TeV 
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Figure 8: The Mjj and Mwjj distributions of Pt,Q>t Wttt — >• and backgrounds at 

the LHC for yfs = 8 TeV and J Cdt = 20 fb -1 . Augmented ATLAS-like cuts as described in 
the text are employed. The open red histograms are the unsealed ttt and pT signals. 



are shown in Fig. [8] for f Cdt = 20 fb -1 . The same cuts as above are used. Counting events 
in the range 120 < Mjj < 170 GeV gives S/y/B = 10.2 for this luminosity but still only 
S/B = 0.050. Despite this large "significance", we remain uncertain of the ability of the 
£ujj channel to settle the questions of CDF's dijet excess and our interpretation of it. 

4. The p^a^ Ztt^ mode 

In view of this situation with the Wttt signal, observation of the isospin partner p^,a^ 
ZtTt of the Wtt^ mode can provide the needed test of the LSTC interpretation of CDF's Wjj 
signal. At the LHC, we predict <r(p^, — >> Zn^) = 2.8 pb, lower than Wn^) = 

4.1 pb because of the reduced phase space, oc p 3 . Then, 

4 Ztt t £ + £~jj) = 190 fb 
for I = e and /i, of which, 80% is due to the p^. This rate is about 10% of the Wttt —> £vdj 
signal. We might expect, therefore, that ~ 10 times the luminosity needed for the Wttt 
signal would be required for the same sensitivity to Zti t . Actually, the situation is better 
than this because there is no QCD multijet background nor $ T resolution to pollute the Zjj 
data. 

Figure [9] shows the Zn T signal and its background, almost entirely from Z + jets, for 
\fs = 8 TeV and J Cdt = 20 fb -1 . The cuts used here are: two electrons or muons of 
opposite charge with p T > 30 GeV, \r) £ \ < 2.5, 80 < M e + e - < 100 GeV and p T (Z) > 50 GeV; 
exactly two jets with p T > 30 GeV and \rjj\ < 2.8; Pr(jj) > 40 GeV, Arj(jj) < 1.75; and 
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Figure 9: The Mjj and M^j distributions of p^, — >■ Z7r^ — >■ £ + £ jj and backgrounds at the 
LHC for y/s = 8 TeV and / Cdt = 20 fb -1 . The cuts used are described in the text. The 
open red histograms are the 7Tt and pr signals. 
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Figure 10: The AR and Ax distributions for p^ — >• Z7r^ — >• £ + £ jj and backgrounds at the 
LHC for yi = 8 TeV and / Cdt = 20 fb -1 . The cuts used are described in the text. The 
open red histograms are the signals. 
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Q = M Z jj — Mjj — M z < 60 GeV. This Q-cnt is very important in reducing the background. 
However, it excludes the 20% of Zttt that comes from production^] These give S/y/B = 
6.2 and S/B = 0.11 for the dijet signal in 120 < M 5j < 170 GeV. The figure also shows the 
M Zjj distribution; it has S/\fB = 6.4 and S/B = 0.12 for 250 < M Zjj < 320 GeV. These 
signal-to-background rates and the position of the dijet signal on the falling backgrounds 
are similar to those in Ref. [2]. Therefore, if our interpretation of the CDF dijet excess is 
correct, both n T —> jj and Pt —> t + £~jj will be observable soon. 

Figure [To| shows the AR and Ax distributions for p T Zn T — >> £ + £~jj. The skyscraper- 
shaped Ax distribution is especially interesting. The background peaks at Ax — 2.3, and 
appears rather symmetrical about this point except that its high side falls more rapidly 
above 2.7 because (Ax)max = tt- The signal's Ax distribution sits atop the skyscraper, 
concentrated in about 330 events in three bins at Ax = 2.2-2.4, whereas the theoretical 
(Ax)min = 2cos _1 (^) = 2.31 for p T —> Zn T . This is just as expected when jet reconstruction 
is taken into account; see Fig. [5} If the actual Ax data, with our cuts, has the shape of our 
simulation, we believe the signal excess can be observed. Similar remarks apply to the shape 



and observability of the slightly broader AR distribution in Fig. 10 



5. The p^a^ WZ mode 

Finally, the decay channel p^ , a? ~~ ^ W ± Z furnishes another important check on the LSTC 
hypothesis provided that sinx ^ 1/4. The dominant contribution, p T —> WlZ l , has an 
angular distribution oc sin 2 9 so that the production is fairly central. We expect g{pt ) clt 
WZ)/a(p T ,a T —> Wn^) ~ (p(Z) / p(tt t )) 3 tan 2 x- The Pythia rates are roughly consistent 
with this. For our input masses and sinx = (§, |, |, |), we obtain the following cross sections: 

a(p T) a T ^WZ trtvi) = (9, 15, 26, 54) fb, (11) 
a(pr, a T ^WZ -> tt'jj) = (27, 48, 80, 170) fb , (12) 
<t(Pt, a T ^WZ -> ivjj) = (90, 155, 260, 555) fb , (13) 
a(p T: a T WW -> tvjj) = (140, 220, 380, 795) fb , (14) 
a(pr, a T Zn T tt'jj) = (205, 200, 190, 145) fb , (15) 

for I = e, p. 

The pt, &t —> t^^t^vi mode has been discussed in Refs. [191 EQ]. It has the advantages 
of cleanliness and freedom from jet uncertainties (except $t resolution). Standard-model 
WZ production at the LHC peaks at M wz = 300 GeV (45], near M PT) and this is the 
dominant background to the 3£is signal. The D0 collaboration searched for this channel 
using the standard LSTC parameters including sinx = 1/3, and excluded it at 95% C.L. up 
to M PT ~ 400 GeV so long as the Pt —> Wttt channel is closed [46] . 



9 We considered Q < 80 GeV to include the a^, but found that the background increased substantially 
faster than the signal. The pt,o>t —> WZ — >> process is included in this simulation, but it also is 

removed by the Q-cut. 
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Figure 11: Left: CMS WZ -> 3£v cross section limits for / Cdt = 4.98 fb" 1 at y/s = 7TeV. 
The LSTC limit curves for sinx = |, |, | assume that = 0.75M PT — 25GeV. Right: 
Two-dimenional exclusion plot for LSTC with sinx = 1/3 as described in the text. The 
CDF mass point is marked by the star. From Ref. [?]. 



The CMS Collaboration recently reported a search for a sequential standard model W 
and for p T , — >• WZ 3£v using 4.98 fb -1 of 7TeV data [?]. The cross section limits and 
M PT vs. Mn T exclusion plot are shown in Fig. [llj The LSTC limit curves for sinx = |, §, \ 
assume that M^ T = 0.75M PT — 25GeV. This stringent assumption significantly enhances 
B(p± -> WZ) above its value for the CDF mass point. For the 2-D exclusion plot, standard 
LSTC parameters, including sinx = 1/3, were used. The CDF mass point is indicated by 
the star. We predicted 21 fb for the signal at 7TeV. Applying a fc-factor of 1.36 in this mass 
range, CMS excludes M nT > 140 GeV at the 95% CL. for M PT = 275 290 GeV. The 95% 
upper limit on the cross section at M PT = 290 GeV is about 20 fb. Using the CMS fc-factor, 
we estimate that the CDF point is allowed for sinx ^ 0.30. 

The dominant background to pr,^r — >• WZ — >> £ + £~jj is Z + jets. As can be inferred 
from Fig. [6] for Wjj production with ATLAS/CDF cuts, the signal will sit at the top of 
the Mjj spectrum. This is what makes the dijet signal in WW/WZ —> ivjj so difficult 
to see. On the plus side, since the LSTC and standard model diboson processes have very 
similar production characteristics, the two signals can be seen with the same cuts and will 
coincide. We simulated this mode and found a promising set of cuts to extract the W jj 
signal. The basic cuts used for the Zjj signal in Sec. 4 were adopted except that we required 
Pr(Z) > 100 GeV, pr(jj) > 70 GeV and 110 < Q = M Zjj - M w - M z < 150 GeV. This 
removed some of the contribution for which the nominal Q = 148 GeV. The mass 
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L = 20 ft) 1 , 8 TeV L = 20 fb" 1 , 8 TeV 




Figure 12: The Mjj and M Z jj distributions of p^ , WZ £ + £~jj and backgrounds at 

the LHC for yfs = 8 TeV and J Cdt = 20 fb -1 . The cuts used are described in the text. The 
open red histograms are the ttt and Pt signals. 



L = 20 fb" 1 , 8 TeV L = 20 fb" 1 , 8 TeV 




Figure 13: The AR and Ax distributions of p^,a^ WZ £ + £~jj and backgrounds at 
the LHC for yfs = 8 TeV and J Cdt = 20 fb -1 . The cuts used are described in the text. The 
open red histograms are the signals. 
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for / Cdt = 20 fb" 



The LSTC signal 
distribution 



distributions for sinx = 1/3 are shown in Fig. 12 

more than doubles the number of standard model W jj events in the Mjj 
and it appears that the dijet signal should be observable with such a data set. Including 
the standard diboson events gives S/y/B = 4.0 and S/B = 0.08 for 60 < M 5j < 100 GeV 



The M Z jj signal is problematic, but it may be possible to combine its significance with that 
for p T Zttt —> £ + £~jj> The AR and Ax distributions are in Fig. 13, The narrow LSTC 
signal and the diboson contribution both peak very near (Ax)min = 2cos~ 1 (vw) = 1-21 and 
they should be observable if the dijet excess is. The £+£~ jj signal is only 60% as large at 
sin x = 1/4 as it is at 1/3. It will be challenging to see it with 20 fb at 8TeV. 



Acknowledgments 

We are grateful to K. Black, T. Bose, P. Catastini, V. Cavaliere, C. Fantasia and M. Mangano 
for valuable conversations and advice. This work was supported by Fermilab operated by 
Fermi Research Alliance, LLC, U.S. Department of Energy Contract DE-AC02-07CH11359 
(EE and AM) and in part by the U.S. Department of Energy under Grant DE-FG02- 
91ER40676 (KL). KL's research was also supported in part by Laboratoire d'Annecy-le-Vieux 
de Physique Theorique (LAPTh) and the CERN Theory Group and he thanks LAPTh and 
CERN for their hospitality. 



19 



Appendix: Nonanalytic Threshold Behavior of da/d(AR) 



1. Kinematics 

We recall first the definition of the angles 0, 0*, 0* and the relevant coordinate systems. 
Choose the z-axis as the direction of the incoming quark in the subprocess cm. frame (or 
the direction of the harder initial-state parton in the pp collision). In the Pt (or a^) rest 
frame, 6 is the polar angle of the ttt velocity i?, the angle it makes with the z-axis. Define 
the xz-plane as the one containing the unit vectors z and i?, so that v — x sin 9 + z cos 0, and 
y — z x x. Define a starred coordinate system in the ttt rest frame by making a rotation 
by angle 9 about the y-axis of the Pt frame. This rotation takes z into z* = v and x into 
x* = xcosO — zsin#. In this frame, let p[ be the unit vector in the direction of the jet 
(parton) making the smaller angle with the direction of v. This angle is 0*; the azimuthal 
angle of p\ = —p* 2 is 0*: 



The jets from i\ T decay are labeled j = 1, 2 and they are assumed massless. Let Ci = + 
and (2 = — ? and c$ = cos0, s$ = sin0, etc. The boosted jets in the lab frame are 



cos 9 — z - cos0* = p* • t), tan0* = p\ y */p\, 



(16) 



|M ffT 7(l + CiUC9.), 
\M^{y + Cjc e »)(a;s e + zq?), 

|M WT Cj((xce - zs g )s0'C^ + ys e *s r ) 



(17) 



_ 1 

where 7 = (1 — ^ 2 ) 2 . 



We want to find the minimum of AR = ^(Arj) 2 + {A(f>) 2 as a function of c e , c e * and c^*. 



From Eq. (17), 




1 + VCp* + (V + Cq*)cq - 7 1 S *C ( f ) *So 
1 + ^* - (V + Q?*)c# + ^Sq+C^Sq 





and 



cos(A0) 



(19) 



^ - ( c i* «i + 7 ( c i 4* + 4*)) ~ 2 7 1 so*c e *s e c c r 



[v 2 s 2 e + (c e *s e + 7 1 s e *c <j) *ce) +(7 1 s e *s <l) *) ] - 4v 2 s 2 e (c *s + 7 1 sq*c <I) *cq) 



1/2 ' 
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2. Minimum of AR 



It clearly is hopeless to deal with the analytic expression of AR as a function of c#, Co*, c^*. 
However, there is a simple way to bypass it. The quantity 

M 2 

A = ^— = coshfAr/) - cos(A0) , (20) 

2PT1PT2 

with Ar] > and < Acj) < tt, is a monotonically increasing function of Ai?. This is seen 
by parametrizing 

Ar] = AR cos A , Acj) = AR sin A (21) 
with A > and A < tt/2 if AR < tt or A < sin-^Tr/Ai?) if Ai? > tt. Then 

5A 

= cos A sinh( Arj) + sin A sin( Acj)) . (22) 

o[AR) 

This is non-negative. It vanishes only for (1) AR = 0, which means A = 0, and this cannot 
happen by its definition, Eq. (20), and for (2) At? = 0, Acj) = tt meaning AR = 7r; the latter 



is a saddle point. This is the "Col du Delta", but it is one-sided, as shown in Fig. 14 




log(l+A) 



Figure 14: The function ln(l + A) defined in Eqs. p0p2| ). The Col du Delta at A = tt/2, 
AR = tt is approached along the road A = tt/2. One cannot go over the pass and down the 
other side for the border is impassable. One must keep climbing along the ridge of increasing 
AR or return via the approach road. 



Minimizing AR thus amounts to minimizing A, which in turn, amounts to maximizing 
PtiPt2- This is much simpler to examine than the original problem. We first maximize 



21 



Vt\Vt2 at fixed c^*, then maximize it with respect to c$*. Since p?j = yPjo ~V% an d Pjo 
depends only on c^*, Pti and ]>T2 are separately maximized at fixed Co* when p\ z = p 2z = 0. 
This requires c$ = sqc^* = 0. Then PtiPt2 = (|7^tt t ) 2 (1 ~~ ^ 2 cos 2 #*) is maximized at 
cq* = 0. In conclusion, AR is minimized if and only if 

Co c e * = c r = . (23) 

This corresponds to two distinct, isolated points in the angular phase space (0* = 7r/2, 37r/2). 
The degeneracy of the minimum is only discrete. At Ai?'s minimum, Arj = and Ac/) = 
cos _1 (2^ 2 — 1) = 2cos _1 (^) = (Ax)min, so that 

( Ai?) min = ( A X ) min = 2 cos" 1 (v) . (24) 
3. Local behavior around cos 6 = cos#* = cos0* = 

We now investigate the behavior of AR as a function of c#, c^* and Co* around its minimum 
at Co = Co* = c^* = by means of a Taylor expansion of at most second order in any of these 
variables. From, Eqs. ( T8fl9 ), we obtain 



(A??) 2 = 4 7 - 2 4, + 0(c 3 ) , (25) 
cos(A0) = cos(Ax) min - (1 - cos(Ax) min ) v 2 (c 2 e + c 2 ,*) + (1 + cos(Ax) min ) 7" 2 c^ + 0(c 3 ) . 

Interpreting the latter equation as: 

cos(A0) = cos(Ax) min - sin(Ax) mi „ (A0 - (Ax) mi „) + 0{{A<f> - (Ax) mi „) 2 ) (26) 

we identify 

A(f> = ( Ax) min + [v 2 tan(( Ax) min /2) (c 2 + c% ) - 7 " 2 cot(( Ax) min /2)c 2 , + 0(c 3 )] . (27) 
Then 

AR = V(A?y) 2 + (A0)2 = (A x ) min + | (6 c 2 + ^c 2 , + &^.cj.) + C(c 3 ) , (28) 

where 

b e = b e * = 2^ 2 tan((A X ) min /2) = 2 V1 ~ l , 

b r = 2 7 - 2 (2/(Ax) mi n-^7)- (29) 

The shape of the surface AR = f \cq , Co* , c^*) in the neighborhood of the minimum 
AR = Axmin is a convex paraboloid with ellipsoidal section whose eigen-directions are 
parallel to the axes of the coordinates c#, Co* and c^*. The curvature is > along each of 
these axes for all < v < 1; i.e. there is no flat direction, as expected from the fact the 
minimum is at isolated point (s). 
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4. Calculation of the singular part of da/d(AR) 

The differential cross section for qq — y — y W/Zn T) followed by ttt —> qq is P*| 



da = 



da(qq -> W/Ztt t ) 
dco 



B(tt t -> qq)) 



dY(Ti T — >► qq) 



r(7r T -> qq) dc e * dc^* 



dco dc$* dcf* . (30) 



To compute the distribution in a compound variable £, such as Ax or Ai?, we use a Fadeev- 
Popov-like trick 

1 = J d(5((-f(c ,c e *,c r )) . (31) 

where / \cq , c$* , c^*) gives the expression of £ in terms of the phase space variables. The 
^-distribution is then 



da 
d( 



= J da(from Eq. (gj) 5 (C - / (c*, c,*c^*)) • (32) 
Let £ = Ai? be slightly above and close to (Ax) m in, and define uj = AR — (Ax)min to shorten 



expressions. Solving Eq. (28) with respect to cq* gives 



Co* 



-l(b e c 2 d + b r c 2 r )+0(<? 



(33) 



Notice that Eq. (33) has to be supplemented by the restriction 



00 



+ Vd+0(c 3 ))>O. 



(34) 



Substituting 

8(^R-f{cg,C0*,C^*)) = (&0.C0*)" 1 [8 (eg* - C e *) + 8 (c e * + Co*)] 6 



u ~ \ ( b <> c d + v4* + °( c ?)) 

(35) 



in Eq. (31) and integrating over c$* leads to the following threshold behavior for the cross 
section: 



f da \ 



threshold 



da(qq -> W/Zn T ) 
dco 



B(n T -> qq) (36) 

c e =c e *c,p* =0 

e[u-l(b e 4 + b r cl,+0(c*))} 



[u-\{b e cl + b r c 2 r +0^))] 



1/2 • 



10 Since there are two points in the (cq, c#*, c^*) phase space where AR has a minimum, = 0*= tt/2 and 
(jf = 7r/2, 37r/2, it is more convenient to use the variable c^* instead of (jf . This introduces (a) the Jacobian 
(1 — c^*) -1 / 2 which is one at c^* = 0; and (b) a factor of two to account for the contributions of the two 
minima in the calculation of the normalization coefficient. 
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It is convenient to trade c#, c^* for new variables p, k\ 
pcosK = ^b /2c , psinft = Jb r /2c r , 



(0<p<^J, 0</^<2tt). (37) 



The integral in Eq. ( |36| ) then yields our final result, the square-root behavior of da/d(AR) 
at threshold: 



da \ 



threshold 



23/2 j AR-(A x )r 
4 6# £u* 



(38) 
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